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Abstract 

We propose a spin model with a new kind of ferromagnetic interaction, which may 
be called ferromagnetic coupling with threshold. In this model the contribute of a 
given spin to the total energy has only two possible values and depends on the number 
of parallel spins among its nearest and next to the nearest neighbors. By mapping 
the model onto the Ising version of the isotropic eight vertex model, we obtain some 
evidence of a low temperature phase made of alternate parallel pluses and minuses 
polymers. 



PACS numbers: 05.50.+q (Ising problems); 64.60.-i (General studies of phase transitions). 
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1. Introduction. 



Ising-like spin systems are deeply studied statistical models which have a surprising richness 
of critical behaviour (see, e. g., [|l|). The axial next nearest neighbor Ising (ANNNI) model 
is the simplest model which describes modulated structures; experimental evidences for such 
phases are furnished by binary alloys (see, e. g., and 0). Modulated ordering in ANNNI 
model is the effect of the competition between ferromagnetic couplings (between nearest 
neighbor spins) and antiferromagnetic couplings (between second-neighbor spins along one 
lattice direction). This competition leads to series of commensurate and incommensurate 
modulated phases of arbitrarily long wavelength (see, e. g., [Q). 

In this paper we show that a system where the spins tend to be aligned to the majority 
of neighbor spins has a low temperature striped phase; this phenomenon is not pointed 
out by the ANNNI model and might lend some insight to the problem of the formation of 
striped-like patterns. 

The model we propose here is an Ising model where the contribution of a given spin to 
the hamiltonian is —(3 if the spin is aligned to the majority of its neighbors and +/3 if it is 
not. In the following we give the definition of the model that we are going to discuss. Let 
us consider the infinite bidimensional square lattice A = and define a Z2 spin variable 
Si G {—1, +1} on each site i E A. We denote hj s E Q = {—1, +1}'^ a configuration of the 
system and define the formal hamiltonian 
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where /3 is a positive real number and 

9 

(j)i{s) = sign Sij} e A, (2) 

where we have denoted by Sij j = 1, ...,9 the nine spins in the 3x3 squared block Bi 
centered on the site i. All the equilibrium properties of the model can be obtained from the 
partition function Z\{P) = Z^sgh exp(— 

We remark that (3 > induces a sort of ferromagnetic coupling among the spin variables, 
but model (0) is ferromagnetic in a different fashion with respect to the Ising model Hj = 
~PJ2<i,j> -Si-Sj, where the sum runs over the pairs of nearest neighbors sites in A. We call 
this kind of coupling ferromagnetic coupling with threshold. 

Indeed, let us suppose Sj = +1: in the Ising model the contribute of the spin Si to 
the total energy of the system decreases when the number of nearest neighbor plus spins 
increases, hence by increasing the number of plus spins among Si nearest neighbors we get 
more and more energetically preferable situations. In our model the contribute of the spin 
Si to the total energy has only two possible values: +/3 if the number of plus spins in Bi is 
less than five, —(3 if it is greater than or equal to five; hence, from the point of view of the 
spin Si, the best situation is reached when four of the eight spins in Bi \ {i} are equal to +1 
and there is no further energetic gain if the number of plus spins in Bi \ {i} increases up to 
eight. 

As it will be explained below, due to the freedom in the choice of number and location 
of the plus (minus) spins around a fixed plus (minus) one, our model shows a very peculiar 
behaviour at low temperatures. In the following we give some evidence of a critical phase 
transition to a low temperature phase, which is a sort of "polymer phase" made of alternate 
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stripes with fluctuating boundaries. 

The occurrence of a striped phase does not depend on the geometry of the blocks Bi, 
which appear in the deflnition of model (|I]). We flnd a striped phase also for the model 

H{s) = -I3 Vsefi , (3) 

where f3 is the inverse temperature and 

5 

4>i{s) = sign {J^Sij} WieA , (4) 
i=i 

where we have denoted by Sij j = 1, ...,5 the flve spins in the "cross" block Ci containing 
the site i and its four nearest neighbors. 

The hamiltonian (|l]) is not strictly ferromagnetic, we mean that it does not satisfy the 
conditions under which Griffiths's inequalities hold 0. Indeed by expanding H as a lin- 
ear combination of products of the spin variables it can be shown that in model (|l|) the 
couplings between nearest and next to the nearest neighbor spins are ferromagnetic, but 
antiferromagnetic terms are present, as well; in the case of model the coupling between 
nearest neighbors is ferromagnetic, while the coupling between next to the nearest neighbors 
is antiferromagnetic. In both the two models a competition between ferromagnetic and an- 
tiferromagnetic couplings occurs, this is similar to what happen in the case of the ANNNI 
model. 

Due to this competition it is not surprising to find striped phases in our model. We 
remark that this competition is essentially an "entropic" effect. The spin Si prefers to have 
just four parallel spins in \ {i}: this is energetically equivalent to the situation with eight 
parallel spins in Bi \ {i}, but the former can be realized in many different ways, while the 
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latter just in one way. The phase transition has been studied by a Renormahzation Group 
Transformation known as "majority rule" (see, e. g., 0,0]); this is suggested by the structure 
of hamiltonian (|l]) . 

The model (|l|) has also an interesting interpretation as a "model of atoms", indeed, it 
can be seen as a model of a binary alloy, namely a mixture of two species of atoms (say A 
and B), such as the Ising model (see, e. g., [§]). But in the present case the configurations 
which are energetically preferred are those in which an A (B) atom has at least four A (B) 
atoms among its nearest and next to the nearest neighbors, no matter if the number of such 
atoms is greater than four. Due to this fact the ordered phase is an alternate sequence of 
stripes of A and B atoms. 

The transformation (^) is widely used for the recovering of noisy images and it is known 
as the median filter (see, e. g., 0). It follows that the hamiltonian ([l|) provides a stochastic 
version of the median filter. 

Finally, we observe that in our model each spin is coupled to a boolean function of its 
neighbors spins; this fact resembles the dynamical rules in the random networks of automata, 
where each spin is updated according to a random boolean function of its neighbors {finite 
dimensional case) or to a random boolean function of a certain number k of random spins 
of the whole network {mean field case, see |llO|). 

Random network of automata exhibit dynamical phase transitions which can be studied 



by the "distance method" (see, e. g., ||T0|)- Our results suggest that, in the finite dimensional 
case, equilibrium phase transitions may arise if one properly chooses the statistical ensamble 
of the boolean functions. 
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The paper is organized as follows: in the next section the expansion of hamiltonians (0) 
and @ in terms of products of spin variables is performed. In section 3 the phase transition 
is studied by mapping the model onto the Ising version of the isotropic 8-vertex model. 
Section 4 summarizes the conclusions. 

2. Some properties of the hamiltonian. 

In this section we discuss some properties of the model introduced above; in particular we 
describe its ground states, we give a rough evaluation of its residual entropy and we expand 
its hamiltonian as sums of products of the spin variables. 

Model (JlD has an infinite number of ground states, which are the configurations satisfying 
the constraint 



such configurations are called median roots, that is configurations invariant under the median 
filtering. 

A complete characterization of bidimensional median roots is still missing, we refer to 
0] for a related discussion and describe them heuristically: a median root is made of two 
typical kinds of structures, namely cluster-like structures and polymer-like structures. In 
Fig. la we show the smallest cluster which is a root, a 12-spins cluster called minroot while 
in Fig. lb a typical polymer structure is depicted. 

It follows that the model has a residual entropy; we obtain a rough estimate of it by 




(5) 
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considering the model defined on a cylinder of 6 x sites with N ^ oo. On a 6 x strip 
(with periodic conditions in the 6-sites direction) we can classify the median roots by the 
configurations on the last two columns. Taking into account the symmetries of the problem 
(parity s — *• — s and invariance under rotations of the cylinder) leads to 30 not-equivalent 
classes. A 30 x 30 transfer matrix Ta^b can be defined as the number of "a" roots on the 
6 X (A^ + 1) cylinder obtained by any "b" root on the 6 x A^ cylinder by adding a column. 

The largest eigenvalue of the transfer matrix provides the residual entropy. We do not 
report here the details, but only quote the result. The largest eigenvalue is found to be 
Aq = 3.57477 and the residual entropy for site |lnAo = 0.21232. To our knowledge there is 
not any estimate in literature for the entropy of the median roots to compare with ours. 

In order to expand the hamiltonian (|l]) as a linear combination of products of the spin 
variables we recall that if f{s) is function of the Z2 spin variables Sj where i ranges over 
some finite set V, then f{s) may be written in a unique way as 

m=j:c{x)s{x) , (6) 

X 

where s{X) = Hiex Si, the sum runs over all the subsets X C V and the numbers c(X) are 
given by 

c{X) = E ^(^)/(^) > (7) 

^ .e{-i,+i}^ 

where \ V\ is the cardinality of V, that is the number of sites in V. 
Now, we observe that H can be written as follows 

H{s) = -f3j:Ms) , (8) 
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where the function /j(s) = Si(j)i{s) is defined on the finite sets Bi Vz G A. Hence, the function 
fi{s) can be expanded as in equation (^: one has to calculate 2^ — 1 = 511 coefficients (that 
is one has to consider all the possible subsets of Bi excepted for the empty set). This number 
could be reduced taking properly into account the symmetries of the hamiltonian. 

We have performed a computer assisted calculation of all the coefficients and, by working 
out the sums in equation (||), we have obtained 

^(^) = -|i:7.E^(c) , (9) 

where the families of clusters Fj and the related coefficients 7j are given in Table 1. We 
remark that each Fj represents a family of a certain kind of cluster of A contained in a 3 x 3 
block B; for example, Fi is the family of all the pairs of nearest neighbor sites, Fg is the 
family of all the plaquettes. 

This calculation shows that model ([^) is not ferromagnetic (see Table 1) in the sense that 
it fulfills the hypothesis under which Griffith's inequalities hold. This model is characterized 
by ferromagnetic and antiferromagnetic couplings: we observe that the 2-spins interactions 
are ferromagnetic if two adjacent columns or rows are involved (nearest and next to the 
nearest spins couplings), while they are antiferromagnetic if the involved columns or rows 
are at distance two lattice spacings. As in the case of ANNNI model, one may expect that 
the result of this competition could be a low temperature striped-phase. 

By performing the same calculation in the case of model (|]) we have obtained 



24 SiSj - 8 Yj ^i^j + 

<i,j> 



«<i,j»> '^i,j,k,i Oij.fc,; 



(10) 



where the five sums run respectively over the pairs of nearest neighbors, the pairs of next to 
the nearest neighbors, the pairs of second neighbors along the lattice directions, the 4-sites 
cluster containing a site i and three of its four nearest neighbors (T-shaped clusters), the 
square 4-sites clusters with sides at 45° with respect to the lattice directions. 

Even model (|^) exhibits a competition between ferromagnetic and antiferromagnetic 
couplings, hence we expect a low temperature striped phase in this case, as well. 

3. The phase transition. 

The structure of the hamiltonian (P suggests to introduce a new set of dynamical variables 
in order to investigate the phase diagram of our model. 

Suppose to partition A into 3x3 squared blocks Ba where a is the site that is in the 
centre of B^, the collection of all these sites a is denoted by A'. On each of these sites a we 
define the new variable 

4 = 0, Va e A'; (11) 

we have defined the new variables by "integrating" over the fiuctuations of the old ones 
on squared 3x3 blocks: we expect that the details of the configurations of the system on 
the scale of three lattice spacings are inessential to describe the peculiarities of the ordered 
phase. 

We have, then, introduced a new model defined on the lattice A', with space of configu- 
rations fi' = { — 1, +1}^'; the equilibrium (unnormalized) measure of this new model is given 
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by 

li\s') = Y.Z{s\s)e~^^^^ , (12) 

sen 

where it has been introduced the probabihty kernel Z{s',s) = DcigA' 5s^,0a(s)- The formal 
hamiltonian of the new model is 

H'{s') = -log fi'{s') + const Vs' G fi'. (13) 

In order to work out the sum in equation (|12D we use the method of the cumulant 
expansion (see [|lT|) writing the hamiltonian (Q) in the form H{s) = Hq{s) + V{s) Vs G fl, 
where Ho{s) = —f!^J2a&A'^a(pa{s) contains the interactions between spins within the same 
block Ba Va e A', while V{s) = ga\A' Si(j)i{s) contains the interactions between spins 

belonging to different blocks. By truncating the cumulant expansion at the first order, one 
can show that the new hamiltonian is in the form 

- H\s') = const + J<> ^'J, + •^«» E « + JuT. «4< , (14) 

<"7> <<«7>> T^Q* 

where the three sums run respectively over all the pairs of nearest neighbors, next to the 
nearest neighbors and over all the plaquettes in A'; we remark that the new model is the 
Ising version of the isotropic eight vertex model With a computer assisted calculation 
we have obtained the new couplings J<>, J<<>> and as functions of the original couphng 
/? (see Fig. 3); their structure is Ja(/9) = /3-Fa(e^^), where Fa are rational functions. 

The variables transformation ( |ri|) and the calculation of the hamiltonian of the new model 
amount to perform the "majority rule" Renormalization Group Transformation (see P, |^). 
We remark that for [3 large enough, with a single step of renormalization, the model (0) is 
mapped into the eight vertex model (|T^ with "negative" nearest and next to the nearest 
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neighbors couplings; that is, starting from a model which is in some sense ferromagnetic, we 
obtain a model with antiferromagnetic couplings. 

The phase diagram of the eight vertex model is well known (see, e. g., ^, |12|): at 
J<<>> < and |J<>| < 2 |J<<>>| this phase diagram is characterized by a critical surface 
separating the paramagnetic and the super antiferromagnetic (SAF) phases, respectively the 
high and the low temperature phases. We observe that as /? increases from zero to infinity 



the renormalized eight vertex model (|T3) undergoes a second order phase transition from 
the paramagnetic to the SAF phase. 

Hence, it there exists a value Pc such that at (3 = (3c our original model is "critical" ; the 
high temperature phase is the usual paramagnetic phase: what about the low temperature 
one? In terms of the new variables s'^ a G A' this phase can be characterized by a staggered 
magnetization m' defined as the difference between the magnetization on even and odd 
columns: m' = in the paramagnetic phase, m' 7^ in the SAF phase. But a column of the 
new model corresponds to a strip of length three lattice spacings in model then, at low 
temperature, model (|I|) exhibits an ordered phase made of an alternate sequence of parallel 
pluses and minuses polymers, infinitely long and large three lattice spacings (see Fig. 2). 
This suggests that the median roots are entropically dominated by polymer-like structures 
developing parallel to one of the lattice directions. 

We have studied model (^) following the scheme used for model (|l|); we have partitioned 
the lattice A as in Fig. 4. The hamiltonian associated to the block variables is 

- H\s') = const + J<> « + J«» E « (15) 

<a7> <<«7>> 

with J<> and J<<>> depending on [3 as in Fig. 5. 
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Even for model @ the renormalized model exhibits a low temperature SAF phase. Due 
to the geometry of the cross blocks, in this case the polymers of the low temperature phase 
are parallel to the dashed hues depicted in Fig. 4. 

4. Conclusions. 

We have proposed an Ising-like spin model with a new kind of coupling in order to point 
out that a simple request, such as the tendency of each spin to be aligned to the majority of 
its neighbors, leads to low temperature striped-phases. 

The expansion of the hamiltonian in terms of products of spin variables pointed out the 
similarities between our model and ANNNI model; this suggests the hypothesis that the 
model undergoes a phase transition to a striped-phase. 

We studied the model by Renormalization Group methods and at the level of approxima- 
tion we have considered, first order of cumulant expansion, the phase transition is observed. 

The parallel polymerization appears to be characteristic of a ferromagnetic coupling with 
threshold. In the spirit of the "atomic" interpretation, this implies that the clusters of A 
atoms, as well as B atoms, turn to be one-dimensional at low temperature, i. e. they loose 
one dimension. 

Our estimations of critical couplings are /3c ~ 25 for model (||) and /3c ~ 5 for model (||). 
The cumulant expansion is known to be relevant near the critical point. A low-temperature 
series expansion would better describe the model in the limit /3 — >■ oo; this is not a trivial 
task, since a rigorous analysis of the statistical properties of the ground states is needed. 
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This will be the topic of further work. 
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Figure Captions 



Figure 1: In Fig. la the smallest cluster invariant to the transformation (§) is shown: black 
squares represent plus spins on a minuses background or, viceversa, minus spins on a pluses 
background. In Fig. lb a portion of a polymer is depicted. The polymer is supposed to be 
infinitely long; it satisfies the constraint (|^), as well. 

Figure 2: The typical low temperature pattern for our model. The average distance between 
polymers of the same sign is six lattice spacings. The constraint (^) is almost satisfied. 

Figure 3: Renormalized couplings as functions of (3 in the case of model ([|). Solid, dashed 
and dotted lines represent respectively J<>, J<<>> and J^. 

Figure 4: The partition of the lattice used to study model (|^). The renormalized variables 
are defined on the blank circles, which form the new lattice. 

Figure 5: Renormalized couplings as functions of [3 in the case of model (|]). Solid and 
dashed lines represent respectively J<> and J<<>>. 
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Table Captions 



Table 1: The coefficients 7^, which are present in equation (|^), are hsted. In the third column 
the families of clusters Fj, to which each coefficient 7j is related, are briefly described. In the 
fourth column a typical cluster C G Fj is depicted: the grid represents the 3x3 block B in 
which ( is contained, the sites of B belonging to C are represented by the black circles. We 
remark that just an example of cluster belonging to Fj is depicted in the fourth column; for 
example, in the case i = 10 one can consider the following clusters, as well: ^3 , ^5 . 
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Table 1 



1 


7* 


QfcibCI ipilOIl OI i j 


xypicai Cluster in i j 


1 


+200 


pairs of nearest neighbors sites 


ffl 


2 


+240 


pairs of next to the nearest neighbors sites 


1 — 1 — 1 

ffl 


3 


-60 


pairs of second-neighbor sites along the lattice 
directions 


ffl 

1 1 1 


4 


-40 


2-sites clusters with sites at distance -\/5 lattice spacings 


ffl 


5 


-20 


pairs of second-neighbor sites along the lattice diagonals 


ffl 


6 


-80 


plaquettes 


ffl 


7 


-40 


4-sites clusters containing the center of the block B and 
occupying a 3 x 2 rectangular block 


ffl 


8 


-20 


4-sites clusters containing the center of the block B and 
occupying the whole B 


ffl 


9 


+24 


4-sites clusters not containing the center of the block B 
and occupying a 3 x 2 rectangular block 


ffl 


10 


+12 


4-sites clusters not containing the center of the block B 

and occupying the whole B 


£a 


11 


+12 


6-sites clusters containing the center of the block B and 
occupying the whole B 


ffl 


12 


+24 


6-sites clusters containing the center of the block B and 
occupying a 3 x 2 rectangular block 


ffl 


13 


-20 


6-sites clusters not containing the center of the block B 
and, necessarily, occupying the whole block B 


ffl 


14 


-20 


8-sites clusters containing the center of the block B 


ffl 


15 


+140 


8-sites clusters not containing the center of the block B 


ffl 
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